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ABSTRACT 

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using 
a sample of mm-selected sources as seen at 250, 350 and 500 \xm by the SPIRE instrument 
on HerscheL From a sample of 63 previously identified mm sources with reliable radio 
identifications in the Great Observatories Origins Deep Survey North and Lockman Hole 
North fields, 46 (73 per cent) are found to have detections in at least one SPIRE band. We 
explore the observed submm/mm colour evolution with redshift, finding that the colours of 
mm sources are adequately described by a modified blackbody with constant optical depth 
r = (p/nwo)^, where ft = -hi. 8 and vq = c/ 100 \im. We find a tight correlation between 
dust temperature and IR luminosity. Using a single model of the dust temperature and IR 
luminosity relation, we derive photometric redshift estimates for the 46 SPIRE-detected mm 
sources. Testing against the 22 sources with known spectroscopic or good quality optical/near- 
IR photometric redshifts, we find submm/mm photometric redshifts offer a redshift accuracy 
of | Az\/(l H- z) = 0.16 ((| Az\) = 0.5 1). Including constraints from the radio-far-IR correlation, 
the accuracy is improved to |Az|/(l + z) — 0.15 ((!Az|) = 0.45). We estimate the redshift 
distribution of mm-selected sources finding a significant excess at z > 3 when compared to 
~850 \im selected samples. 

Key words: galaxies: high-redshift - galaxies: statistics - submillimetre: galaxies. 


1 INTRODUCTION 

As is the case in the local Universe , star formation at high redshift 
takes place in regions rich in molecular gas and dust. Franceschini 
et al. (1991) and Blain & Longair (1993) thus anticipated that the 
distant galaxies hosting the majority of this activity would be dis- 
covered in the submm/mm waveband, which benefits from a strong, 
negative -correction for the rest-frame far-infrared (far-lR) por- 
tion of their spectral energy distributions (SEDs). Observing at 450 
and 850 p.m, the Submillimetre Common-User Bolometer Array 
(SCUBA) submm camera on the James Clerk Maxwell Telescope 
(Holland et al. 1999) made it possible to detect the first examples 
of this galaxy population and represented a fundamental turning 
point in our understanding of galaxy formation in the distant Uni- 
verse (Snail, Ivison & Blain 1997; Barger et al. 1998; Hughes et al. 
1998; Dey et al. 1999; Eales et al. 1999). 

Early work on the submm galaxy population was often hindered 
by small numbers of sources, and significant difficulties in find- 
ing multiwavelength counterparts given the relatively low angular 
resolution (full width at half-maximum, FWHM > lOarcsec) of 
ground-tased submm/mm facilities. This situation has been much 
improved by a large increase in the number of submm/mm facil- 
ities, larger scale submm surveys [e.g. SHADES, Mortier et al. 
2005; LAB OCA Extended Chandra Deep Field South Submm Sur- 
vey (LESS), WeiB et al. 2009], improvements in the quality of 
overlapping multiwavelength data (in particular radio and mid-IR 
imaging. , as well as advances in cross -matching algorithms (Ivison 
et al. 2007; Chapin et al. 2009, 2011; Roseboom et al. 2009). As a 


result, our understanding of these galaxies has accelerated in recent 
years (e.g. Alexander et al. 2005; Chapman et al. 2005; Greve et al. 
2005; Pope et al. 2008; Swinbank et al. 2008; Tacconi et al. 2008; 
Men6ndez-Delmestre et al. 2009). 

The advent of the SPIRE instrument (Griffin et al. 2010) onboard 
Hcrschel 1 (Pilbratt et al. 2010) promises to shed new light on many 
issues relating to submm galaxies. In photometry mode, SPIRE 
operates in three submm bands centred on 250, 350 and 500 p.m, 
and thus the large-scale surveys conducted with SPIRE will detect 
many 10000 s of submm galaxies (e.g. Eales et al. 2010; Oliver 
et al., in preparation ). SPIRE observations of existing submm/mm 
survey fields provide the crucial SED information shortward of 
500 \im which cannot be routinely accessed from the ground (cf. 
Kovdcs et al. 2006), allowing us to characterize their IR luminosity 
(Li R ) and dust temperature (7p) with unprecedented accuracy (e.g. 
Chapman et al. 2010; Magnelli et al. 2010). 

Armed with this new information, we can re-evaluate the poten- 
tial of using submm (and mm) wavelength data alone to estimate 
the redshifts of submm galaxies. As submm galaxies tend to be 
both high redshift « z ) ~ 2; Chapman et al. 2005) and highly dust 
obscured, they are extremely faint at optical/near-IR wavelengths 
and hence obtaining redshifts via optical spectroscopy is challeng- 
ing. Because of this, the prospects for submm/mm photometric 


1 Herschel is a European Space Agency space observatory with science 
instruments provided by Principal Investigator consortia. It is open for pro- 
posals for observing time from the worldwide astronomical community. 
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redshifts have been discussed many times since the early SCUBA 
surveys (Blain 1 999; Hughes et al. 2002; Aretxaga et al. 2003; Blain, 
Barnard & Chapman 2003; Pope & Chary 2010; Schulz et al. 2010), 
with mixed results. The main obstacle for submm/mm photometric 
redshifts is that redshift and dust temperature are degenerate for the 
simplest, single temperature, modified blackbody SED. To break 
this degeneracy, some other constraint is needed, such as the IR 
luminosity to dust temperature relation (Blain 1999; Blain et al. 
2003) or information at other wavelengths, physically related to the 
far-IR emission (e.g. radio or mid-IR; Carilli & Yun 1 999; Aretxaga, 
Hughes Sc Dunlop 2005; Aretxaga et al. 2007). 

In this paper, we investigate the potential of submm/mm photo- 
metric redshifts using a sample selected at SPIRE and mm wave- 
lengths. We focus on two key survey fields: Great Observatories 
Origins Deep Survey North (GOODS-N) and Lockman Hole North, 
which contain some of the deepest imaging on the sky at mm, 
submm and radio wavelengths. In particular, the Owen & Morrison 
(2008) and Morrison et al. (201 0) Very Large Array (VLA) imaging 
of the Lockman Hole and GOODS-N fields represent the deepest 
‘blank-field’ images currently available at 1 .4 GHz. The depth of 
the radio imaging is crucial as without the positional information 
from identifications in high angular resolution radio imaging it is 
often not possible to correctly deblend the highly confused SPIRE 
imaging. 

In Section 2, we introduce the parent mm- wavelength catalogues 
utilized in this work. Section 3 describes the process undertaken to 
find mult: wavelength identifications for the mm and SPIRE sources. 
Section 4 presents the SPIRE detection statistics for our mm- 
wavelength sample, and describes the properties of those sources 
with both mm and SPIRE detections. Finally, Section 5 considers 
the scientific implications of these results, and we present our con- 
clusions in Section 6. Throughout we assume a A cold dark matter 
cosmology with = 0.7, Q ra = 0.3 and Hq = 70 km $ -1 Mpc -1 . 

2 DATA 

The starting point for this study is the available mm-wavelength 
imaging in the GOODS-N and Lockman Hole North fields. 

In GOODS-N, we make use of the combined MAMBO-AzTEC 
image and source catalogue of Penner et al. (2011). Combining 
the pre-existing 1.1-mm AzTEC (Perera et al. 2008) and MAMBO 
1.2-mm (Greve et al. 2008) imaging of GOODS-N, the resulting 
image has a typical noise of 0.5 mJy at an effective wavelength of 
1.16 mm over an area of 0.08 deg 2 . This combined image yields 
a list of 41 1.16-mm sources detected at a significance of >3.8 <t 
and 5u6nim = 2-10 mJy. In Lockman Hole North, we make use 
of the recent MAMBO 1.2-mm imaging and source catalogue of 
Lindner et al. (2011). Lindner et al. present a list of 41 1.2-mm 
sources detected with signal-to-noise ratio (S/N) > 4 and Si. 2 mm = 
2-5 mJy. 

In addition, we make use of the SPIRE imaging at 250, 350 and 
500 pm obtained as part of the Science Demonstration Phase of Her- 
schelby the Herschel Multi-tiered Extragalactic Survey (HerMBS; 2 
Oliver et al., in preparation). The SPIRE instrument, its in-orbit per- 
formance and its scientific capabilities are described by Griffin et al. 
(2010); its calibration methods and accuracy are outlined in S win- 
yard et al. (2010). We use images produced as described in Levenson 
et al. (2010). The SPIRE instrument has an FWHM of 18, 25 and 
36 arcsec at 250, 350 and 500 fim, respectively. HerMES imaging 

2 hermes.sussex.ac.uk 


of GOODS-N and Lockman Hole North has a typical 5er depth at 
all SPIRE wavelengths of —5 and ~10mJy, respectively, ignoring 
the dominant contribution from confusion (Smith et al. 2012). 

At radio wavelengths, we make use of the VLA 1.4-GHz source 
catalogue from Morrison et al. (2010) and Owen & Morrison 
(2008) for the GOODS-N and Lockman Hole fields, respectively. 
In GOODS-N, the Morrison et al. (2010) imaging has a sensitivity 
of 3.9 pJy beam -1 with the sensitivity declining to 8 pJybeam -1 
at a radius of 15arcmin from the centre. In Lockman Hole 
North, the Owen & Morrison (2008) imaging has a sensitivity of 
2.7 pJy beam -1 over the central 40 x 40 arcmin 2 region overlapping 
the mm and SPIRE observations. 

3 MULTIWAVELENGTH IDENTIFICATIONS 
FOR SUBMM/MM- DETECT ED SOURCES 

The large beam size (> 10 arcsec FWHM) of typical submm/mm 
imaging facilities means that accurately compiling multi wavelength 
identifications and photometry for submm/mm-detected sources is 
often a significant challenge (e.g. Lilly et al. 1999; Ivison et al. 2007; 
Roseboom et al. 2009; Chapin et al. 2011). This challenge arises 
mainly in two ways; first, the large beam size, and low S/N, of typical 
submm/mm surveys means that catalogued sources will have quite 
large positional uncertainties. This makes matching to catalogues at 
other wavelengths difficult, as there may be more than one potential 
counterpart within the positional errors of the submm/mm source. 
Indeed, for faint submm sources at high redshift, the source density 
of matching catalogues deep enough to contain the true match may 
be so high that more than one match can always be found. 

The second problem is source confusion; more than one astro- 
nomical source is present in the submm/mm beam, and hence cata- 
logued sources are made up of multiple astronomical objects. This 
problem is particularly worrying as it affects the integrity of the 
submm/mm source catalogues themselves. 

Several statistical methods which can account for one or both of 
these issues exist (Downes et al. 1986; Roseboom et al. 2009, 2010; 
Chapin et al. 2011). For the data sets considered here, one or both 
of these issues are present. Source identification is a major issue for 
the mm-wavelength data sets, but source confusion is not, as mm 
sources have large positional uncertainties, but low source densities. 
The SPIRE data sets have both large positional uncertainties and 
high source densities, and so are equally affected by both. 

Thus, different approaches are needed to reliably produce mul- 
tiwavelength associations for the different submm/mm data sets. 
In practice, this means that it is generally not possible to reliably 
identify SPIRE counterparts to mm-detected sources without first 
making identifications in some other catalogue with higher spatial 
resolution (i.e. 1.4 GHz or 24 vim). To emphasize this point, the 
typical flux density of our mm-detected sample, after correcting 
for flux boosting, is S i. 2m m ~ 2 mJy. Taking rough estimates of the 
mean dust temperature and redshift of submm/mm sources, 35 K 
(Kov&cs et al. 2006; Chapman et al. 2010) and z = 2.5 (Wardlow 
et al. 201 1 ), and assuming a single temperature modified blackbody 
SED, gives predicted SPIRE fluxes of 13, 19 and 15 mJy at 250, 350 
and 500 vim, respectively. These flux levels are below the confusion 
‘limit’ at these wavelengths, with the best estimate of the noise from 
confused sources cr C0nf = 6 mJy (Glenn et al. 2010; Nguyen et al. 
2010). Hence, recovering the SPIRE counterparts to mm-detected 
sources will require a technique which can reduce the effect of con- 
fusion, by using prior information at other wavelengths (Roseboom 
et al. 2010; Chapin et al. 201 1). Putting the issue of source confu- 
sion aside, even if these sources could be recovered reliably, source 
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identification would still pose a serious challenge. The positional 
uncertainties of both the mm- and SPIRE-detected sources (apo* ~ 
2-3 arcsec) require the use of a 9-12 arc sec matching radius in or- 
der to recover the bulk of the true mm-SPIRE associations. Given 
the kno v /n number density of SPIRE sources at these flux levels 
(~7000 leg -2 ; Glenn et al. 2010), we would expect between 14 and 
25 per cent of these associations to be chance alignments. 

In this work, we first match the mm-detected source lists to 
1.4-GHz sources, using the well-established catalogue based cross- 
matching techniques of Downes et al. (1986) and then use these 
1.4 GHz source positions, along with the source positions of all 
known 1.4 GIIz and 24 pm sources in these fields, to ‘deblend’ 
the SPIRE maps into their individual contributors and thus provide 
the best estimate of the SPIRE photometry at the locations of our 
mm-selected sample. 


3.1 1.4 GHz identifications of mm-selected sources 

In the GOODS-N field, we search for potential 1 .4 GHz counter- 
parts in the Morrison et al. (2010) VLA source catalogue to the 41 
1.16-mm sources from Penner et al. (201 1) using a lOarcsec search 
radius. 28 of the 41 sources have at least one potential counter- 
part within this search radius. To determine the reliability of these 
matches, we estimate the probability of a chance alignment using the 
P-statistic, as defined by Downes et al. (1986). Defining a reliable 
match tc be the one with a less than 5 per cent probability of being 
a chance alignment, 32 reliable matches are found for 24 objects, 
with eight sources having more than one reliable radio counterpart. 

For the sources with more than one reliable counterpart, we 
take the counterpart with the lower P-statistic to be the correct 
identification. 

One difficult case for our matching is HDF 850. 1 (ID 14 in the cat- 
alogue of Penner et al. 2011), the brightest submm source detected in 
the original 850 pm SCUBA imaging of GOODS-N (Hughes et al. 
1998). Our naive cross-matching approach would associate this ob- 
ject with a 56 pjy radio source located 5.3 arcsec away. However, 
subsequent near-IR and high-resolution Submillimeter Array ob- 
servations show that the true position is much closer to the original 
SCUBA position (Dunlop et al. 2004; Cowie et al. 2009), although 
no distinct 1.4 GHz radio source can be found near this position in 
the Morrison et al. (2010) image. Here we adopt the Dunlop et al. 
(2004) position for matching with the SPIRE data. 

Other difficult cases are those sources for which the submm 
emission is later shown, via high-resolution mm-wavelength inter- 
ferometry, to originate from multiple components. In these cases, 
the 1 .4 GHz identification(s) is often found to lie coincident with one 
(or more) of these components (e.g. Wang et al. 2011). Amongst 
our sample are two examples of this, GN20 (G1 in our sample; 
Daddi et al. 2009a; Carilli et al. 2010, 2011) and GN21 (G19 in 
our sample; Wang et al. 2011). In both cases, we have simply taken 
a single radio identification, the one with the lowest P-statistic, as 
being solely responsible for the mm flux density. For Gl, the ra- 
dio identification corresponds to GN20 from Daddi et al. (2009a), 
while for G19 this corresponds to GOODS 850- 13c in the list of 
Wang et al. (2011). While it may turn out that a large fraction of 
mm-detected sources do indeed contain multiple components, in 
both Gl and G 19 the single radio identification we have chosen 
is found to be responsible for >50 of the submm/mm flux (Daddi 
et al. 20C9a; Wang et al. 201 1). Importantly, our SPIRE photometry 
pipeline accounts for all 1.4 GHz, mm and 24 pm detected sources, 
so it should not be affected by this issue. 


Encouragingly, for the 28 sources from Penner et al. (201 1 ) which 
were previously identified in the AzTEC and/or MAMBO surveys, 
we recover the same identifications made by both Chapin et al. 
(2009, hereafter C09), for the 25 sources in common, and Greve 
et al. (2008), for the eight sources in common. 

In the Lockmann Hole North field, we make use of the existing 
matches between the 1 .2-mm sources and 1.4 GHz radio catalogue 
presented in Lindner et al. (201 1). 40 of the 41 1.2-mm sources 
are found to have a 1 .4 GHz counterpart within 8 arcsec, with 39 
of these found to have less than a 5 per cent probability of being a 
chance alignment and hence are deemed reliable. 

It is worth noting that three of these reliable radio counterparts 
are not found in the published Owen & Morrison (2008) 1.4 GHz 
source lists, and come from a deeper extraction performed in the 
vicinity of the MAMBO sources by Lindner et al. (2011). For these 
sources, we adopt the quoted 1.4 GHz positions and flux densities 
from Lindner et al. (201 1). 

Thus, we are left with a sample of 63 mm-selected sources with 
high angular resolution 1.4 GHz positions with which to match to 
the SPIRE data. 

Interestingly, there is a significant disparity between the identi- 
fication rates in the two fields considered here. Both Penner et al. 
(2011) and Lindner et al. (2011) list 41 mm-detected sources with 
a similar noise^ la ~ 0.7 mJy, although in both cases the noise 
varies significantly across the field. However, the effective area of 
the two surveys is quite different. The Penner etal. (201 1) combined 
AzTEC/MAMBO image of GOODS-N covers an area of 0.08 deg 2 , 
while the Lindner et al. (2011) MAMBO image of Lockman North 
covers an area of 0.16 deg 2 . Thus, it is clear that for the Penner 
et al. GOODS-N sample to find the same number of objects similar 
to one-half the area, it must be identifying a fainter population of 
sources. The differential number counts at these wavelengths are 
known to be well described by a power law <W/d5ocS*, with an 
exponent <5 ~ —3. So to observe a difference of a factor of 2 in 
surface density, the GOODS-N sample must be effectively a factor 
of V2 deeper, ignoring the ^-correction between 1.16 and 1.2 mm. 

3.2 SPIRE photometric method 

The most difficult obstacle in determining accurate SPIRE pho- 
tometry is the effect of source confusion, i.e. the contributions of 
numerous faint sources within a single SPIRE resolution element, 
centred on the target of interest. For this reason, several authors have 
developed techniques that utilize the positions of sources detected 
at other wavelengths, usually 24 pm and 1 A GHz, to disentangle 
the various contributions from discrete sources to the SPIRE flux 
in a given beam element (e.g. Roseboom et al. 2010; Chapin et al. 
2011). This process is made possible by the high correspondence 
between the 24-pm and 1.4-GHz populations and those observed 
at far-IR wavelengths; >80 per cent of the cosmic IR background 
at SPIRE wavelengths can be accounted for by 24- pm sources with 
S 2 4 > 25 pJy (Marsden et al. 2009; Pascale et al. 2009), while the 
strong correlation between the far-IR and radio luminosity is known 
to hold across a wide range in redshift and luminosity (e.g. Ibar et al. 
2008; Ivison et al. 2010). 

Here we follow the prescription presented by Roseboom et al. 
(2010, hereafter R10) to measure the SPIRE fluxes of mm sources 
in GOODS-N and Lockman Hole North, with a number of modifi- 
cations. Below, we summarize the R10 methodology, then highlight 
our departures from it. 

R10 assume that the SPIRE map can be fitted directly by assuming 
that the positions of all sources contributing significantly to the 
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map are known (i.e. previously detected at 24 pm, 1.4 GHz or other 
wavelengths), and that only the SPIRE flux density of each of these 
sources is unknown. Algebraically this means the observed SPIRE 
map, d, can be described as a series of point sources with fluxes,/, 
at known positions, x. Thus, the observed map can be described as 

d=P x f + & , 

where P x is the point response function (PRF) centred at posi- 
tion x and 8 is some unknown noise term. As discussed in RIO, 
this equation can be easily inverted, using linear algebra, to ob- 
tain the best-fitting values for the fluxes. However, one potential 
pitfall of this approach is that of overfitting. If many of our in- 
put sources are intrinsically faint and the source density is high, 
this approach will give poor results. RIO solve this issue by using 
model-selection techniques to eliminate those sources that are un- 
necessary to describe the map. Specifically, the Akaike information 
criterion 'Akaike 1974) is used to filter the input list of sources via 
backward elimination. This is very time consuming and hence the 
number of sources that could be fitted simultaneously was limited 
to -100-200 in RIO. 

Here, we replace this model-selection approach with the adap- 
tive Least Absolute Shrinkage and Selection Operator (LASSO; 
Tibshirani 1996; Zou 2006). In particular, we combine a non- 
negative application of this technique - the Non-Negative Least 
Squares path (NNLSpath) algorithm of ter Braak et al. (2010) - 
with the adaptive weighting scheme suggested by Zou (2006). In 
Zou (2006), the weights are tunable, based on initial estimates of 
the parameters; here we use the probability of a chance alignment, 
<t > , as defined in R10. Specifically, <j> is calculated separately for 
the three samples which make up our input list: radio, 24 p.m and 
mm-detected sources. In cases where a source is detected in more 
than one of these bands, the lowest value of (p is adopted. As mm- 
detected sources have the lowest areal density, these typically have 
higher weightings, followed by radio sources and then 24 pm. 

Accuracy in model selection is the primary advantage. The adap- 
tive LASSO is known to have 'oracle' properties 3 with the correct 
choice of weightings; however, there is also a significant advantage 
in computational speed. With this approach, we can solve for large 
numbers of candidate sources (up to 10000) in a realistic amount 
of time. 

In addition, we modify the background estimation. Rather than 
using an unconstrained local estimate of the background, the back- 
ground is estimated globally, assuming a flat background across the 
total extent of the map. This approach is preferred as it is more 
robust against incompleteness in the prior input list, and because 
there are no significant large-scale correlations in the noise present 
in our SPIRE maps due to the exceptionally low-frequency ‘knee’ 
to the 1// noise (Griffin et al. 2010). 

Finally, potential new SPIRE sources, i.e. those that were previ- 
ously unknown at 24 p.m or radio wavelengths, are considered by 
examining the residual maps after fitting the SPIRE fluxes of known 
sources. 

3.3 Estimating the true noise for SPIRE photometry method 

Accurately estimating the effect of confusion on an individual 
source is non-trivial. As the SPIRE maps have a mean intensity 

3 In model selection, a method is said to have oracle properties if it can 
reliably identify which free parameters are truly present in the underlying 
system, and which are unnecessary or spurious. 


of zero, for blind source catalogues we can expect that the typical 
confusion noise is given by the width of the distribution of flux den- 
sities [P(Z))] in a PRF-convolved image. For SPIRE, we know these 
values are 5.8, 6.3 and 6.8 mjy beam' 1 at 250, 350 and 500 pm, 
respectively (Nguyen et al. 2010). However, our source extraction 
process attempts to find the true zero-point of the map, and simul- 
taneously fits the positions of sources found at other wavelengths, 
in an attempt to reduce this noise. However, it is difficult to as- 
sess for a given source how much we have ‘resolved’ the confusing 
background, and how much is left as unresolved fluctuations in the 
map. 

One simple way to assess the global improvement found by 
our method is to use completely synthetic maps and catalogues 
produced using input mock catalogues known to give reasonable 
agreement to the known far-IR and submm number counts (e.g. 
Femandez-Conde et al. 2008). Testing against simulated SPIRE 
maps with the same properties as our GOODS-N data (i.e. instru- 
mental noise, number counts and pixel intensity distribution), we 
find a typical total noise, i.e. confusion and instrumental noise, of 
cr t 0 tai = 2.8, 3.9 and 3.7 mjy at 250, 350 and 500 pm, respectively. 
These values compare quite favourably to the SPIRE photometry 
process described in R10, which gave <r to tai = 4.9, 5.4 and 6.5 mJy 
at 250, 350 and 500 pm for the same simulations. 

While testing the performance of our method on simulated data 
sets is of interest, we really need an accurate estimate of the true 
error on the SPIRE photometry of our actual sources. One approach 
would be to estimate the global confusion noise in our SPIRE pho- 
tometry by investigating residual intensity fluctuations in the SPIRE 
maps, after subtracting the contributions from ‘resolved’ sources. 
This approach is appealing, as we know that the residual maps con- 
tain information about sources not included in our prior source list, 
as well as any artefacts from incorrect astrometry and/or assump- 
tions about the PRF. 

However, there is one complication for this approach; variations 
in sensitivity in the prior source lists across the SPIRE map will 
mean the ability to resolve confusing sources varies spatially. Thus, 
we must assess the confusion noise on scales small enough that the 
prior source list sensitivity is invariant. This is particularly difficult 
in this work, as the combination of 24 pm and 1.4 GHz prior source 
catalogues results in large variations across the field. This is due to 
both the difference in relative sensitivity of 24 pm and 1.4 GHz to 
SPIRE sources, as well as the tapering of the 1.4 GHz sensitivity 
towards the edges of its coverage. 

We utilize a variant of this approach to estimate the true error 
for our confused SPIRE sources here. Specifically, we construct 
PRF-convolved residual maps for each of the SPIRE bands, using 
our best estimate of the SPIRE photometry at the point of each 
input source and the known SPIRE beam. Then, at the position of 
each source we are interested in, we assess the total noise per pixel 
by measuring the standard deviation of surrounding pixels within 
a 20-pixel radius (2, 2.8 and 4arcmin for the SPIRE 250, 350 and 
500 pm maps used here). 

A further complication arises in converting this noise per pixel 
into the noise for a point source as we know that the confusing 
background must be correlated, on small scales by the PRF and on 
larger scales by the intrinsic clustering of the sources. This means 
that in calculating the total noise estimate for a given point source, 
we must take into account the uncorrelated instrumental noise in 
each pixel and the correlated noise from confusion. For simplicity, 
we assume that the correlated noise on small scales is dominated by 
the PRF, and assume that the shape of the PRF gives the covariance 
between pixels directly, i.e. if we assume the SPIRE photometry is 
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estimated by 

«. E idjPhl 

LW 

where d\ is the intensity in pixel /, P, the PRF for the point source 
in pixel i and a, the instrumental noise in pixel i, the noise estimate, 
including contributions from both confusion noise and instrumental 
noise, can be given by 

o&rt = (P T C-'P)-', 

where C is the covariance matrix, the diagonal elements of which are 
given by 1/cr- 2 , and the off diagonal elements i,j given by Pj/Ci/aj, 
where F is the PRF assuming a point source centred at position i. 

Using this approach, we find the typical la total noise for our 
mm-detected sources is = 2.2, 2.6 and 2.5 mJy in the central 
part of the GOODS-N field for the 250, 350 and 500 pm bands, re- 
spectively, increasing to 3.7, 4.8 and 5.7 mJy ala radius of 20 arcmin. 
Similarly, in the Lockman North field, the total noise is found to be 
3.5, 3.6 and 4.7 mJy at the centre of the field, and 4.4, 4.7, 6 mJy at 
a radius of 20 arcmin. 

It is worth noting that these variations are solely due to the vari- 
ations in the depth of the prior source list; the SPIRE sensitivity is 
uniform across both fields within the coverage of the mm observa- 
tions (Levenson et al. 2010). 


4 RESULTS 

4.1 SPIRE detection statistics 

We define the SPIRE measurement as a detection if it has >3<r 
significance in the flux density estimate and the goodness of the 
photometric fit to the map is high (reduced x 2 < 5, in a 15 x 15 pixel 
region centred on the source position). The x 2 cut is necessary to 
remove cases where bright or unaccounted for sources, or residual 
non-Gaussian noise (e.g. Galactic cirrus or uncorrected thermal 
drifts) introduce signal that cannot be accounted for by our input 
source model. 

Table 1 presents our SPIRE photometry for all mm-selected sam- 
ple in GOODS-N, while Table 2 presents SPIRE photometry for 
mm sources in the Lockman Hole North field. 

In GOODS-N, 17 of the 24 1.16-mm sources with good identifi- 
cations are detected in at least one SPIRE band, while in Lockman 
North 30 of the 39 sources with identifications are detected in at 
least one SPIRE band. These figures give a marginally lower de- 
tection rate in GOODS-N (71 per cent) as compared to Lockman 
North (77 per cent), although this is consistent with the slightly 
deeper nature of the GOODS-N mm imaging (see Section 3.1). 

As GOODS-N and Lockman Hole North are well-established 
survey fields, a large number of spectroscopic, and good quality 
photometric, redshifts are available. In GOODS-N, nine sources 
have spectroscopic redshifts, either from targeted observations of 
submm galaxies and radio sources by Chapman et al. (2005, and in 


Table 1. Measured SPIRE photometry for mm sources in GOODS-N. 3a upper limits are given for sources which are not detected above this significance in 
the SPIRE imaging. For the mm photometry, deboosted flux densities are quoted. The x 2 represents the reduced x 2 in a 15 x 15 pixel region centred on the 
source position. ID number refers to the equivalent value from Pennsr et al. (201 1). C09 refers to the ID numbers of sources in common with Chapin et al. 
(2009). G08 refers to the ID numbers of sources in common with Greve et al. (2008). 1 .4 GHz flux density is from the catalogue of Morrison et al. (2010). 


Name 

C09" 

GOS* 

RA 

(°) 

Dec. 

n 

5250 

(mJy) 

*250 

5350 

(mJy) 

xho 

_ ^ 
11 

*500 

■Sf.2 mm 
(mJy) 

5l.4GHz 

(mJy) 

Gl* 

1 

1 

189.29950 

62.37013 

18.6 ±2.3 

0.97 

31.7 ±2.7 

2.16 

24.3 ± 2.7 

1.10 

10.5 ± 0.7 

0.079 ± 0.014 

G2 d 

3 

2 

189.13918 

62.23579 

18.1 ±2.3 

1.77 

15.2 ± 2.7 

2.11 

<7.8 

1.68 

5.2 ± 0.6 

0.036 ± 0.004 

G3 d 

5 

4 

189.37811 

62.21633 

46.3 ± 2.0 

1.25 

46.5 ± 2.3 

1.75 

32.7 ± 2.3 

1.27 

4.4 ± 0.6 

0.123 ± 0.004 

G4 d 

7 

3 

189.29727 

62.22540 

38.6 ±2.1 

1.76 

33.0 ± 2.4 

2.62 

23.2 ± 2.3 

1.30 

4.0 ±0.5 

0.131 ±0.009 

G5 

2 

13 

189.13301 

62.28749 

<6.9 

1.58 

<8.7 

1.68 

<7.8 

0.88 

4.1 ± 0.6 

0.026 ±0.005 

Gff 1 

6 

5 

189.11346 

62.10158 

9.2 ± 2.4 

1.23 

13.5 ± 3.0 

1.61 

5.2 ± 2.9 

1.14 

4.5 ± 0.7 

0.053 ± 0.012 

C7 

4 

12 

188.95979 

62.17839 

<7.5 

1.53 

<8.1 

2.31 

<8.4 

2.09 

4.6 ± 0.7 

0.028 ± 0.005 

GS d 

26 

6 

189.30772 

62.30739 

8.2 ±2.2 

1.91 

8.1 ±2.8 

3.76 

<7.8 

2.94 

3.1 ± 0.6 

0.652 ± 0.004 

G9 d 

11 

14 

189.14937 

62.11888 

33.9 ± 2.1 

1.45 

20.3 ± 2.7 

1.16 

<8.4 

1.43 

2.9 ± 0.6 

0.046 ± 0.008 

Gll d 

13 

15 

188.97157 

62.22708 

21.0 ±2.4 

1.25 

<9.3 

4.06 

<8.4 

2.39 

2.6 ± 0.6 

0.045 ± 0.005 

G13 

12 


189.13604 

62.10586 

<6.9 

1.44 

<8.7 

1.35 

<8.7 

1.25 

2.7 ± 0.7 

0.027 ± 0.005 

G14 

14 

25 

189.21696 

62.20714 

< 6.3 

1.05 

<7.5 

1.77 

<7.2 

1.31 

2.4 ± 0.6 

0.016 ± 0.004* 

G16 



189.20117 

62.35201 

<9 

1.00 

<9.3 

1.98 

<8.1 

1.24 

2.5 ± 0.6 

0.030 ± 0.005 

G18 rf 

16 


189.06708 

62.25388 

<7.2 

1.35 

<9.3 

2.95 

14.1 ±2.6 

1.46 

2.1 ± 0.6 

0.036 ± 0.005 

G19" 


29 

189.30020 

62.20341 

11.2 ± 21 

1.50 

11.3 ±2.4 

1.28 

17.3 ± 2.4 

1.09 

2.1 ±0.6 

0.032 ±0.004 

G21 d 

9 


189.40941 

62.29355 

7.1 ±2.0 

1.10 

12.1 ± 2.6 

1.53 

<8.4 

1.01 

2.1 ± 0.6 

0.026 ± 0.006 

G24 

24 


189.03583 

62.24325 

<7.5 

1.34 

<9.6 

3.08 

<8.1 

2.36 

2.0 ±0.6 

0.039 ± 0.004 

G26 d 

15 


188.94971 

62.25811 

<8.1 

2.33 

<10.2 

4.22 

10.1 ±2.9 

1.93 

2.1 ± 0.7 

0.124 ±0.006 

G27 d 

18 


189.42151 

62.20591 

31.0 ±2.3 

1.19 

27.9 ± 2.4 

0.99 

30.3 ± 2.5 

0.88 

1.9 ±0.6 

0.033 ± 0.005 

G29 d 

25 


189.21550 

62.08417 

12.6 ±2.9 

1.78 

15.1 ±3.3 

2.57 

17.7 ± 3.0 

1.43 

2.6 ± 0.9 

0.079 ± 0.006 

G30 d 

17 


188.92167 

62.24403 

11.7 ±2.6 

1.59 

<10.2 

3.80 

<9 

2.27 

2.2 ± 0.7 

0.029 ± 0.005 

G31 d 



189.32254 

62.13453 

35.1 ±2.5 

3.20 

32.1 ±3.0 

1.95 

22.5 ± 2.7 

0.96 

2.0 ± 0.7 

0.050 ± 0.008 

G34 



189.14537 

62.32322 

<7.2 

1.36 

<87 

2.28 

<10.3 

0.99 

1.8 ± 0.6 

0.085 ± 0.005 

G41 d 



189.55217 

62.24859 

25.3 ± 3.2 

2.48 

25.6 ± 3.5 

9.72 

15.3 ±3.7 

2.97 

1.9 ± 0.7 

0.056 ± 0.006 


"Typically written with prefix AzGN. 

^Typically written with prefix GN1200. 
c Deboosted flux density as described in Penner et al. (201 1). 
^Detected in at least one SPIRE band. 

*1.4 GHz flux density from Dunlop et al. (2004). 
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Thble 2. Measured SPIRE photometry for mm sources in Lockman Hole North. 3 o upper limits are given for sources which are not detected 
above this significance in the SPIRE imaging. For the mm photometry, deboosted flux densities are quoted. The x 2 represents the reduced x 2 in 
a 15 x 15 pixel region centred on the source position. ID number refers to the equivalent value from Lindner et al. (2011). 1.4 GHz flux density 
is from the catalogue of Owen & Morrison (2008), as presented by Lindner et al. (201 1). 


ID 

RA 

Dec. 

*$250 

xfso 

*$350 

*350 

•$500 

*500 

•$i.2mm (deboosted) 

•$1.4GHz 


O 

o 

(mjy) 


(mJy) 


(mJy) 


(mjy) 

(mJy) 

LL J 

161.75083 

59.01878 

80.1 ± 3.4 

0.92 

63.2 ±3.8 

0.85 

40.6 ± 4.4 

0.89 

3-5^6 

0.279 ±0.015 

L2 r 

161.61192 

59.09578 

<11.7 

1.07 

12 5 ± 4.0 

0.99 

<13.8 

1.08 

3-8ioj 

0.036 ± 0.006 

L3 

161.63112 

58.84889 

<10.8 

0.84 

<11.1 

0.92 

<13.2 

0.87 

3-8^7 

0.459 ± 0.025 

L4 

161.53050 

58.90389 

<10.8 

0.83 

<11.4 

1.03 

<14.1 

0.79 

7 7+0.5 

Z * -0.5 

0.030 ± 0.004 

L5 G 

161.85592 

59.06014 

32.8 ± 3.9 

1.21 

44.5 ±4.1 

1.08 

35.6 ±4.6 

0.92 

4 0 

4U -0.9 

0.051 ±0.006 

L6 n 

161.66112 

58.93681 

52.4 ± 3.6 

0.91 

51.5 ±3.6 

0.75 

17.5 ±4.4 

2.17 

7 7+0.4- 
Z •*- 0.4 

0.160 ±0.010 

L7~ 

161.75054 

58.91150 

21.0 ±3.6 

0.84 

18.6 ±3.5 

0.97 

31.7 ±4.4 

0.78 

7 3+°* 
ZJ -0.5 

0.042 ± 0.005 

L8 n 

161.63779 

58.86642 

30.0 ±3.6 

0.99 

29.8 ± 3.7 

0.94 

16.5 ±4.4 

1.25 

7 <7+0.6 
Z/ -0.6 

0.098 ± 0.010 

L9 

161.77042 

58.83556 

<10.8 

0.96 

<11.1 

0.91 

<13.8 

1.21 

1 Q+1-0 

0.023 ±0.005 

L10 fl 

161.59604 

58.99300 

13.3 ± 3.5 

0.72 

<14.1 

0.94 

19.8 ± 4.4 

0.91 

2.4^j 

0.078 ± 0.007 

Lll rt 

161.48721 

58.88856 

58.2 ± 3.7 

1.02 

56.6 ± 3.8 

1.26 

18.8 ±4.6 

1.44 

2.7^ 

0.315 ±0.019 

L12" 

161.19829 

59.00997 

25.3 ± 3.7 

0.84 

36.0 ± 3.9 

1.11 

24.4 ± 4.7 

0.78 

7 7+0.6 

Z ’ —0.7 

0.273 ± 0.019 

L13 a 

161.53646 

58.97458 

12.1 ± 3.6 

0.98 

17.4 ± 3.6 

0.80 

<14.1 

0.82 

2 l+°- 5 
Z * ] -0.5 

0.198 ±0.010 

L14" 

161.64937 

59.13014 

39.0 ± 3.9 

0.99 

46.2 ± 4.0 

1.07 

33.9 ± 4.6 

1.65 

3 O^* 9 

0.097 ± 0.006 

L15” 

161.86654 

58.87058 

55.8 ± 3.9 

1.01 

69.1 ± 3.6 

0.95 

49.7 ± 4.8 

0.80 

3 0+ 0 * 9 

0.180 ±0.018 

L16" 

161.83633 

58.86472 

<11.7 

1.04 

11.2 ±3.6 

0.90 

<14.1 

0.95 

7 7+0.8 
Z *'-0.8 

0.048 ± 0.007 

L17« 

161.54250 

59.04508 

10.4 ±3.7 

1.15 

14.6 ± 3.9 

1.26 

<13.5 

0.84 

2 0 +OS 
Z ’ -0.5 

0.027 ± 0.003 

LlS fl 

161.73021 

58.83444 

39.9 ± 3.6 

0.81 

37.8 ± 3.7 

0.91 

25.4 ± 4.6 

0.95 

2 

z ’*-l.l 

0.105 ±0.008 

L19 

161.25833 

59.06758 

<11.4 

0.91 

<12.3 

1.23 

<14.7 

1.45 

9 n+0.6 
zU -0.6 

0.029 ± 0.004 

L21" 

161.37575 

59.11008 

<11.7 

1.02 

12.0 ±4 

1.11 

16.0 ± 4.9 

1.03 

1 R+ 0 - 5 
15 -0.5 

0.036 ± 0.004 

L22" 

161.51763 

59.08039 

32.7 ± 3.6 

1.06 

23.7 ± 3.9 

1.01 

<14.1 

1.00 

2 0 +0 6 
z ‘ u -0.6 

0.166 ±0.010 

L23 

161.67129 

58.89050 

<11.1 

0.97 

<10.8 

0.92 

<13.2 

0.72 


0.032 ± 0.005 

L24 a 

161.25279 

59.12600 

11.0 ±3.6 

0.92 

<12.0 

0.76 

17.6 ±4.6 

0.79 


0.068 ±0.009 

L25 

161.41679 

59.06333 

<12.0 

1.25 

<12.6 

1.66 

<15.0 

6.59 

1 7+0.6 

0.024 ± 0.007 

L27 a 

161.76000 

58.85086 

17.4 ±3.5 

0.91 

19.5 ±3.6 

0.71 

<13.8 

1.16 

2 0+ 08 
z * -0.9 

0.078 ± 0.014 

L28 

161.58708 

58.90944 

<10.5 

0.74 

<11.1 

1.25 

<13.5 

1.04 

1 5+ 0 - 5 
1J -0.5 

0.025 ± 0.006 

L29" 

161.48129 

59.15453 

140.5 ± 3.9 

2.47 

64.4 ±3.9 

0.83 

29.1 ±4.8 

1.46 

i o+0.8 
J °-0.8 

0.307 ± 0.043 

L30* 

161.29329 

59.06864 

<11.7 

1.02 

133 ±4.1 

1.28 

<15.0 

1.59 


0.074 ±0.007 

L31 fl 

161.60379 

58.89644 

<10.5 

1.04 

20.8 ± 3.7 

1.08 

13.1 ±4.3 

0.89 

1 -5-2:5 

0.043 ± 0.004 

L32* 

161.41579 

58.90692 

18.6 ±3.8 

1.07 

21.7 ±3.9 

0.98 

17.6 ± 52 

1.08 


0.047 ±0.004 

L33° 

161.39608 

58.84714 

28.3 ± 3.7 

1.19 

32.5 ± 3.8 

0.96 

26.5 ± 5.1 

1,13 

2-7 

0.071 ±0.010 

L34 fl 

161.22346 

58.97819 

<10.8 

0.94 

13.3 ±3.8 

0.91 

<14.1 

0.74 

l-5t? 

0.041 ± 0.008 

L35 fl 

161.82546 

58.92383 

19.7 ±3.7 

0.95 

24.4 ±3,7 

0.97 

18.2 ±4.5 

1,14 

1 5 +0 - 7 

0.060 ± 0.006 

L36 

161.53375 

59.12889 

<11.7 

1.40 

<11.7 

0.98 

<14.1 

1.13 


0.016 ±0.003 

L37« 

161.54575 

58.87914 

52.0 ± 3.6 

1.01 

39.9 ± 3.8 

1.09 

16.8 ±4.5 

0.95 

i c+0.8 

! * 3 -0.7 

0.160 ±0.011 

L38 

161.18704 

59.13836 

<12.0 

0.86 

<11.7 

0.81 

<14.4 

4.57 

I 5 +0 9 
* • -0.9 

0.068 ±0.016 

L39° 

161.55012 

59.04261 

<14.1 

1.12 

<11.4 

1.22 

23.3 ±4.5 

0.81 

1 4 +07 
l * 4 -0.6 

0.034 ± 0.003 

L40“ 

161.74304 

59.10931 

20.4 ± 3.7 

1.09 

12.2 ±4.0 

0.89 

<14.1 

1.05 

1 5 +I - 2 

0.057 ± 0.005 

L41° 

161.50129 

58.91794 

12.1 ± 3.7 

0.86 

11.6 ±3.8 

0.87 

<14.1 

1.19 

1 4 +07 

0.072 ±0.012 


a Detected in at least one SPIRE band. 


preparation) or from the many magnitude-limited redshift surveys 
of the field which have been conducted (see Barger, Cowie & 

Wang 2008, for a recent compilation of spectroscopic redshifts in 
GOODS-NT). In the Lockman Hole North, spectroscopic redshifts 
are available for four sources, mostly from mid-IR spectroscopy of 
24 pm selected sources (Fiolet et al. 2010) or WIYN spectroscopy 
of radio sources (Owen & Morrison 2009). Photometric redshifts 
are available for an additional nine mm sources in Lockman Hole 
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North from the catalogue of Strazzullo et al. (2010). Here we only 
consider photometric redshifts which have a quality flag of ‘A A’, i.e. 
the estimated photo-z error is less than 0.2(1 + z ), from Strazzullo 
et al. (2010). In total, good quality redshifts are available for 22 of 
our parent sample of 63 mm sources. These redshifts are quoted 
in Tables 1 and 2. All of the sources with known spectroscopic 
and/or reliable optical/near-IR photometric redshifts are detected 
by SPIRE. This is not particularly surprising, as those sources with 
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redshift information from the optical/near-IR wavelength range tend 
to be more intrinsically luminous and/or at slightly lower redshift 
(Wardlcw et al. 2011). 

To estimate the reliability of our detections, we test random posi- 
tions in the SPIRE maps. Test positions are required to be >3 arcsec 
from other sources in our input list, and within the region covered 
by our 24- pm and 1 .4-GHz catalogue. This is the typical minimum 
separation of sources in our input list, due to the resolutions at 24 pm 
and 1 .4 GHz. If we impose the same selection criteria on the fluxes 
recovered for these random positions, 2.7 per cent are detected in 
any SPIRE band at a significance of >3 o. 

Taking this as a guide, if our parent sample of 63 mm sources 
were completely spurious, then we would expect to recover SPIRE 
detections in any band for two sources. 


4.2 Subrnm-mm colours 

For the jse of submm/mm photometry, either alone or in combina- 
tion with mid-IR or radio data, to be useful as a redshift estimator, 
there must be a reasonably strong, and unique, matching between 
submm/mm colour evolution and redshift. 

In Fig. 1, we compare' the SPIRE-mm colours to the expected 
evolution for a modified blackbody parametrized by a single dust 
temperamre, J D , and optical depth r = (v/v 0 A i.e. 

S,«{| -«,[-(!)' 

where is the Planck function for a single dust temperature T D . In 
order to give a better representation to the mid-IR SED, we replace 
the blackbody SED on the Wien side of the SED with a power 
law, S v (X v a , where a is fixed to —2 (Blain et al: 2003; Conley 
et al. 201 1). To ensure a smooth SED, the power law is used at fre- 
quencies where dlogS/dlog v < a. As we have a limited number of 
photometric bands to constrain this model, we fix the dust opacity 
and optical depth to values of p = +1.8 and t> 0 = 100 H-m, respec- 
tively. We adopt this definition of a modified blackbody throughout 
this paper. 

The sources with known spectroscopic or reliable optical/near-IR 
photometric redshifts from Tables 3 and 4 are coloured coded by 
redshift. 

While the errors on the SPIRE and mm photometry are large, 
the modified blackbody model can sufficiently explain the observed 
submm-mm colours of the majority of SPIRE-mm sources. This 
is qualitatively in agreement with recent results from Kov£cs et al. 
(2010), although they warn that more sophisticated models (in par- 
ticular, one which takes account of the range of dust temperatures 
likely to occur within a star-forming galaxy) are needed to match 
shorter wavelength data (i.e. less than rest frame 30 pm). A modi- 
fied blackbody redshift and temperature are completely degenerate 
at z < 4; thus, changing the assumed dust temperature will change 
the predicted submm-mm colour as a function of redshift, but it 
will not change the position of the model track in colour-colour 
space. Thus, if the dust temperature were known, the redshift could 
be inferred directly, or vice versa. 

It can be seen that the highest redshift objects are found to have the 
‘reddest* submm/mm colours, i.e. large Shuh/Sspire or S 500 /S 350 flux 
ratios, supporting claims that identifying red, or ‘cold*, SPIRE/mm 
sources is an effective way of isolating very high redshift submm 
galaxies (e.g. Pope & Chary 2010). 
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Figure 1. Submm-mm colour evolution of mm-selected sources: 
S500/S350 versus S 35 o/S 25 o(top), Si.2mm/S350 versus S350/S230 (middle) and 
£i.2mm/£500 versus S 500 /S 350 (bottom). Only sources with >3<r detections in 
each of the relevant SPIRE bands are shown in each panel. The track taken 
by a modified blackbody, with T& = 40 K and optical depth r = (v/voA 
where /J = +1.8 and vo = c/100 pm, is also shown. Sources located in 
the GOODS-N field are shown in squares, those in Lockman Hoie North 
are shown as triangles. Symbols are colour coded by redshift, with filled 
symbols designating those sources with spectroscopic redshifts and open 
symbols those with reliable [Sz < 0.2(1 + z)] optical/near-IR photometric 
redshifts. Black open symbols represent sources with no redshift estimate. 
The submm-mm colours are reasonably well described by our simple mod- 
ified blackbody model. 

4.3 IR luminosity to dust temperature relation for SPIRE-mm 
sources 

As a significant fraction of SPIRE-mm sources have good redshift 
information (22 of 46 SPIRE-detected mm sources), we can inves- 
tigate the IR luminosity to dust temperature relation directly for our 
sample of SPIRE-mm sources. 

For each source with a spectroscopic or optical/near-IR photo- 
metric redshift (excluding the very low z source L29), we fit a 
modified blackbody to the SPIRE and MAMBO photometry (and 
AzTEC and SCUBA where available), assuming a dust emissivity 
of /3 = +1.8, the redshift and the functional form for a modified 
blackbody given in Section 4.2. In the Lockman Hole North field. 
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Table 3. Radio-far-IR photometric redshifts, far-IR luminosities and 
dust temperatures for mm-detected sources in the GOODS-N field. 
Spectroscopic redshifts from the literature are quoted where available. 
Far-IR luminosities (8-1000 pm) and dust temperatures are presented 
fur sources with good redshift information. 


GOODS-N 




ID 

Redshift 

Radio/far-IR photo-z 

I/IrCEq) 

Td( K) 

Gl" 

4.05* 

5.00^5° c 

13.2 ±0.1 

39 ±1 

G2" 

4.05“' 

4.20t“;®5' 

13.1 ±0.2 

43 ±2 

G3" 

2.98' 

2-63*$' 

13.2 ±0.1 

50 ±2 

G4" 

1.996^ 

2 | -34-0.00 c 
Z ' U -0.17 

12.7 ±0.2 

36 ±2 

G5 


4.5* 



G 6 " 


4.07*$' 



G7 


4.6* 



G 8 “ 


o 7 - 1 0.4 h 
-*• 7 “ 0.2 



G9 

0.95 1 

0 qi + 013 c 
uyi -0.15 



Gil" 

2.098' 

a 07 + 0.22 c 
U - yZ -0.41 

12.8 ±1.3 

73 ±36 

G13 


3.95* 



G14 


4.1 ±0.5' 



G16 


3.7* 



G18" 

2.5& 

2 56*9- 18 " 
Z ~ H, -4U7 

12.5 ± 0.6 

35 ±8 

G19" 

2.91/ 

7 Q 1 +A-24 c 

2.91^23 

12.6 ± 0.3 

36 ± 3 

G21" 

3.1?/ 

7 07 +O .34 c 
Z.52_o27 

12.6 ± 0.4 

38 ±4 

G24 


3.1* 



G26° 


O 24+ 0 - 05 c 
az4 -0.09 



G27" 


7 4-5 +0.04 c 
Z ‘ 4:5 -0.18 



G29" 


0 ft 1 +0-26 c 
Z ‘ 8 , -0.3l 



G30" 


9 41 +0-42 c 
Z * 41 -0.26 



G31" 


UZ 3 - 0 . 1 I 



G34 


2 . 2 * 



G41" 


t 07 +O.I 7 C 
1 - 0.21 




"Detected in at least one SPIRE band. 

^Spectroscopic redshift from Daddi et al. (2009a). 

"Redshift estimate from simultaneous fitting of Lir-7d and Lir- 

f'MGHz relations. 

^Spectroscopic redshift from Daddi et al. (2009b). 

" Spectroscopic redshift from Chapman et al. (in preparation). 

f Spectroscopic redshift from Chapman et al. (2005). 

* Redshift estimate from mm-1 .4 GHz spectral slope. 

Redshift estimate from fitting of Lir-7d as radio flux boosted by 

AGN activity. 

'Spectroscopic redshift from Barger et al. (2008). 

-'Radio-far-IR photo-z estimate from Dunlop et al. (2004). 

PACS 100 and/or 160 pm photometry is available from the HerMES 
survey for six sources and is included in the SED fits. In GOODS- 
N, SCUBA 850 pm photometry is available for five sources from 
the catalogue of Pope et al. (2005) and is included in the SED fits. 
IR luminosities are calculated for each source by integrating the 
modified blackbody fit in the range 8-1000 pm. IR luminosities 
and dust temperatures for sources with good redshift information in 
GOODS-N and Lockman Hole North are presented in Tables 3 and 
4, respectively. 

Fig. 2 shows the best-fitting modified blackbody IR luminosity 
versus dust temperature for our sample of SPIRE-mm sources. One 
source, Gil, is omitted as our modified blackbody model is not 
well constrained by the observed data (see Table 3). Also shown 
is the local relationship for IRAS 1.2 Jy sources from Chapman 


Table 4. Radio-far-IR photometric redshifts, far-IR luminosities and 
dust temperatures for mm-detected sources in the Lockman Hole 
North field. Spectroscopic or reliable [&z < 0.2(1 + z)] oplical/near-IR 
photometric redshifts from the literature are quoted where available. 
Far-IR luminosities (8-1000 pm) and dust temperatures are presented 
for sources with good redshift information. 


Lockman Hole North 


ID 

Redshift 

Radio/far-IR photo-z 

Lir(Lq) 

7d(K) 

Ll fl 

2.562* 

7 Ofi +0 09 c 
ZUO -0.18 

13.3 ±0.2 

54 ± 3 

L2" 


4 OQ+ 0 - 27 c 
4 - uy -0.28 



L3 


1.8" 



L4 


4.4" 



L5" 

3 ± 0.42" 

a an+0.07 c 
^ u -0.29 

13.1 ±0.2 

43 ±2 

L6" 

2.03^ 

1 0 1 +0.06 c 

17 

12.9 ±0.2 

43 ±2 

L7" 


2 . 74 *$' 



L8" 


2.07ig;g' 



L9 


3.9" 



L10" 


J.U1_ 0 .21 



Lll" 

1.95^ 

1 71 +O.I 8 c 
L 71 -0.13 

12.9 ±0.2 

41 ±3 

LI 2" 

2.2 ± 0 . 2 " 

9 £ 4 + 0.08 c 
zo*»_0.23 

12.7 ± 0.2 

33 ±2 

LI 3" 


A 07 +O.I 8 c 
UJZ - 0 . 2 I 



LI 4" 

2.3 ± 0.5" 

7 01 +0.16 c 
z - 61 -0.18 

12.9 ± 0.2 

37 ±2 

L15" 

2.8 ± 0 . 2 " 

•* AO -0.16 

13.3 ± 0.2 

45 ±2 

L16" 


3 94 + 0 - 5 J c 
3 Z 4 -0.33 



LI 7" 


2 25 c 
ZO ->_0.19 



LI 8 " 


7 1 c+ 0.11 c 
Z * i:> -0.24 



L19 


4.1* 



L21" 


2.521“^' 



L22" 

1.4 ±0.4" 

ft 7 ft+^ H c 
Ui8 -0.28 

12.3 ± 0.5 

32 ±4 

L23 


3 . 6 " 



L24" 

3.2 ±0.4" 

7 C 7 +O .43 c 
^-3+_0.20 

12.7 ± 0.5 

40 ± 5 

L25 


3.5* 



L27" 


1 62 ri).45 r 
1 ' OZ -0.27 



L28 


3 . 8 " 



L29" 

0.004* 

ft 7 ft +0.28 c 
u, * 5 O -0J3 

8.0 ±0.3 

22 ±2 

L30" 


A 7 1 + 0.20 c 
u - 71 -0.41 



L31" 

2.9 ±0.3" 

2 44 +° 21 c 
Z44 -0.24 

12.7 ± 0.4 

40 ± 4 

L32" 

1.3 ± 0 . 2 " 

7 7 1 +0.25 c 
Z,Z1 — 0.27 

12.0 ±0.4 

24 ±2 

L33" 


/> ^o+0.16 c 
Z - W -0.34 



L34" 


7 4£+0.50 c 
z/K) -0.32 



L35" 


7 14 +O.I 8 c 
Z * 14 -0.24 



L36 


4 . 5 " 



L37" 

1.7 ± 0.2" 

1 . 17 ^:!*' 



L38 


3 . 6 " 



L39" 


ZJ 5 y -o. 2 i 



L40" 


ft 7R+0 - 13 c 
U/JS -0.46 



L41" 


1 4Q+ I-®* c 
1 ,4y — 0.87 




"Detected in at least one SPIRE band. 
b Spectroscopic redshift from Polletta et al. (2006). 
c Redshift estimate from simultaneous fitting of Lir-7d and L \ R - 
Li. 4 GHz relations. 

^Redshift estimate from mm-radio spectral slope as presented by 
Lindner et al.(201 1). 

"Photometric redshift from Strazzullo et al. (2010). 

I Spectroscopic redshift from Fiolet et al. (2010). 

^Spectroscopic redshift from Owen & Morrison (2009). 


© 2011 The Authors, MNRAS 419, 2758-2773 
Monthly Notices of the Royal Astronomical Society © 201 1 RAS 



SPIRE-mm photometric redshifts 2767 



Figure 2. IR luminosity versus dust temperature for those SPIRE-mm 
sources with good redshift information. Source G 1 1 is omitted as a modified 
blackbody model is not well constrained by the observed data. Symbols are 
colour ceded by redshift, with GOODS-N sources denoted by squares and 
Lockman Hole North sources denoted by triangles. The solid line is the local 
Z'IR-T’d relation from the 1 .2 Jy IRAS sample as determined by Chapman 
et al. (2CD3), where we have converted from Sjoo/Sgo to Td assuming a 
modified blackbody with optical depth r = (u/voA with p = +1.8 and 
vq = c/lGO jam. The dashed line shows a modification to this relation given 
by fitting equation (1) directly to our SPIRE-mm sources. The best-fitting 
parameters are y = 0.193, <5 = —0.03 and C* = 0.3. L* is held fixed at 
5 x 10 10 . 


et al. (2003), after converting from IRAS iSW^eo colour to To via 
our definition of a modified blackbody. It is clear that the SPIRE- 
mm sources form a tight sequence in L^-Td space. However, we 
can see that our sample is systematically biased towards colder 
dust temperatures when compared to the local Lir-T d relationship 
(Chapman et al. 2003). This bias towards colder dust temperatures 
for submm/mm-selected samples is now well established (Chapman 
et al. 2010; Magdis et al. 2010; Magnelli et al. 2010). 

The bias towards colder dust temperatures causes the Lir-T d rela- 
tion to be much tighter than for an unbiased sample. To parametrize 
the observed Lir-7d relation for mm-selected sources, we fit the 
functional form given by Chapman et al. (2003) to our mm-selected 
sample. Specifically, Chapman et al. (2003) describe the mean IRAS 
colour via 

logC 0 = logC. + Slog(l + ^) +Kl0g ( , + I^)’ (1) 

where C* = 0.45, S = —0.02, = 5 x 10 10 Lq and y = 0.16, 

while Z/hr is the integrated IR luminosity from 3-1 100 p.m and Co is 
the ratio of the 60-to- 1 00 \im IRAS flux density. We fit equation ( 1 ) to 
our observed L\r-T d relation of our SPIRE-mm sources, adopting 
a conversion between IRAS 60-to- 100 pm flux ratio and T D given 
by a modified blackbody with = +1.8. We leave L* fixed at 
the Chapman et al. (2003) value as we do not have enough low- 
luminosity sources to constrain this parameter. Our best fit is given 
by C* = —0.3, S = —0.03 and y = 0.193. The track taken by our 
best-fitting Li R -7 d relation is shown in Fig. 2. 
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5 DISCUSSION 

5.1 Submm photometric redshifts 

As shown in Section 4.2, there is a reasonable correspondence be- 
tween submm/mm colour and redshift, which is likely a by-product 
of our SPIRE-mm sources demonstrating a strong relationship be- 
tween dust temperature and luminosity, as seen in Fig. 2. Thus, 
if we can assume a relationship between dust temperature and IR 
luminosity, we can break the temperature-redshift degeneracy and 
produce photometric redshifts which rely on submm/mm data alone. 

We explore the potential of such an approach here by making 
use of the 22 SPIRE-mm sources with good redshift information as 
a test set. We assume that the Lir-7d relationship for SPIRE-mm 
sources is given by our best fit of equation (1) as seen in Fig. 2. 
Additional photometry from PACS and SCUBA is not considered in 
the photometric redshift estimation as these are only available for a 
small subset of the sample; we want a fair assessment of the ability 
of SPIRE and mm-wavelength data alone to estimate the redshift. 

Fig. 4 compares our photometric redshift estimates to the known 
redshifts for these sources. The mean photometric redshift accuracy 
is found to be | Az|/(1 + z) = 0.16 with the mean offset ((|Az|) = 
0.51). All of our SPIRE-mm sources are detected in the radio at 

1.4 GHz, so we can hope to improve on our photometric redshift 
estimates by also requiring that the sources lie on the radio-to-far-IR 
correlation (van der Kruit 1971). As our SPIRE photometry is close 
to the peak of the far-IR emission at z ~ 1-3, it is not sufficient to 
assume a spectral index between the SPIRE/mm bands and 1 .4 GHz 
as is typically done for submm-radio photo-z (e.g. Carilli & Yun o 
1999). Here we assume the ratio of Lir to Lmghz is given by a con- 
slant value <7ir = log 10 [(L 1R /3.75 x 10 12 W)/(L 140Hz /W Hz' 1 )] = 

2.4 (Ivison et al. 2010). We also assume that the radio SED is 
adequately described by a power law, S„ocv“, with a = —0.75 
(Ibar et al. 2010). Fig. 3 shows q as a function of redshift for 
the 22 SPIRE-mm sources with good redshifts. Good agreement 
is seen between our sample and the median (and standard devia- 
tion) of q for SPIRE sources as measured by Ivison et al. (2010). 



Figure 3. Radio-to-far-IR correlation as a function of redshift ay 
probed by the qm parameter, where qjR — log io [(Z-ir/3 .75 x 
IQ 12 W)/(Z,i.4ghz/WHz -1 )]. The dashed line and shaded region represent 
the median and ±1 <t values for SPIRE sources from Ivison et al. (2010). 
Squares and triangles represent sources with reliable redshift information 
in GOODS-N and Lockman Hole North, respectively. Good agreement is 
seen with the Ivison et al. (2010) values. No significant evolution with 
redshift is observed, suggesting our single #ir model should introduce no 
redshift-dependent bias to the photometric redshift estimates. 
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Importantly, no evolution with redshift is observed, similar to Ivi- 
son et al. (2010), and hence our single q model should introduce no 
redshift-dependent bias. As none of the sources in Fig. 3 lies sig- 
nificantly off the radio-to-far-IR correlation, we are confident that 
potential contributions to the radio flux density from active galactic 
nucleus (AGN) activity may be ignored. 

To get the best estimate of the redshift, we simultaneously fit 
the L ir-7d and Lir-Li. 4 gh 2 relations, assuming the far-IR emission 
is described by a modified blackbody. After including the radio- 
to-far-IR correlation, the photometric redshift accuracy improves 
to |Az|/(l + z) = 0.15 (or (|Az|> = 0.45). If we consider only 
those sources with spectroscopic redshifts, the accuracy improves 
again to | Az|/(1 -f z) = 0.12 (or (| Az|) = 0.35). Interestingly, this 
same improvement is not seen for the L [R -T D only method when 
the sample is restricted to only those with spectroscopic redshifts 
(I Az|/(1 + z) = 0.16 or ((| Az|> = 0.55)). 

While finding the redshift that gives the maximum likelihood 
redshift (i.e. the best fit) is the primary aim, it is also of interest 
to examine the probability distribution to see how wide a range in 
redshift also gives a reasonable fit. For each source, we calculate the 
redshift range in which 68 per cent of the probability is enclosed. 
These points are shown as error bars in Fig. 4. Using this approach, 
the typical estimated redshift error for photometric redshifts using 
only the L\r-T d is a z = 0.6, while including the radio-to-far-IR 
correlation reduces this to a z = 0.4. Encouragingly these numbers 
are quite similar to the true error measured by comparing to known 
redshifts. 

Many previous works have attempted to use a combination of 
radio and ground-based submm/mm data to estimate photometric 


redshifts. Aretxaga et al. (2007) estimated submm-radio photo- 
metric redshifts for 120 850- pm sources detected in the SHADES 
survey (Mortier et al. 2005), 58 of which had previously known 
spectroscopic or robust optical/near-IR photometric redshifts. They 
make use of two distinct photo-z techniques: a simple 850 fim to 
1.4 GHz spectral index estimator, similar to that proposed by Carilli 
& Yun (1999), and full fitting of the far-IR to radio SED to a set 
of well-defined templates. Comparing the results of these two tech- 
niques to the known redshifts gives a mean photometric redshift 
error of | Az| = 0.9 for the spectral index method and | Az| = 0.6 
for the template fitting method. 

More recently, Biggs et al. (2011) used a similar spectral index 
approach to estimate photometric redshifts for 870- pm sources de- 
tected in the LESS survey of ECDFS. Comparing the 48 sources 
which have spectroscopic redshifts, or reliable 17-band photometric 
redshifts (Wardlow et al. 2011), the photometric redshift accuracy 
is found to be | Az|/(1 + z) = 0.25 or a mean deviation of (| Az|) = 
0 . 68 . 

Thus, it is clear that photometric redshifts estimated using the 
Lir -7 d relation alone offer a slight improvement over radio-submm 
redshifts, while using a combination of the Lir-T d relation and 
radio-far-IR correlation, the accuracy of submm/mm-radio photo- 
metric redshifts can be increased by almost a factor of 2. By con- 
trast, a similar study by Dannerbauer et al. (2010) using PACS 100 
and 160 p.m, in combination with mid-IR, radio and submm/mm 
(but not SPIRE), data find that the accuracy of photometric red- 
shifts using IR templates is significantly worse if PACS data are 
included. This highlights both the need for SPIRE data, which is 
at or close to the peak wavelength of the far-IR emission at these 



Figure 4. Spectroscopic versus photometric redshift estimates for submm/mm photometry alone, assuming a universal Lir-7d relation (left), and combined 
with radio constraints, assuming a universal radio-to-far-IR correlation (right). Filled symbols designate sources with spectroscopic redshifts, while open 
symbols designate those sources with optical/near-IR photometric redshifts only. For submm/mm photometric redshifts, we achieve a mean photo-z accuracy 
of | Az|/(1 4- z) = 0.16 ((| Az|) = 0.51), while if radio constraints are included this is reduced to | Az|/(1 + z) = 0.15 «| Az|) = 0.45). Error bars are calculated 
by finding the redshift range which contains 68 per cent of the probability. 
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redshifts, as well as the need for improved SED models in the mid-IR 
(5-30 pm). 

5.2 Redshift distribution of submin/mm sources 

The redshift distribution of sources selected at different submm and 
mm wavelengths has been a subject of much debate. Many previ- 
ous works have claimed that sources selected at mm wavelengths 
are preferentially found at higher redshifts than those at —850 pm 
(Dannerbauer et al. 2004; Younger et al. 2007; Greve et al. 2008; 
Coppin et al. 2009). In contrast, others have claimed that the 
—850 pm and mm-selected populations are at similar redshifts, with 
the apparent differences originating from selection effects (Wardlow 
et al. 2011). 

Thus, it is of interest to re-examine this using our sample of 
SPIRE-mm sources. While only 22 of the 46 SPIRE-detected mm 
sources have reliable spectroscopic or photometric redshifts, as we 
have seen in Section 5.1, it is possible to estimate the redshifts of 
SPIRE-mm sources with reasonable accuracy using a combination 
of subrrni/mm and radio data alone. 

As the submm photometric redshift errors are large, and our sam- 
ple of SPIRE-mm sources small, we prefer to use the full probability 
distribution function (PDF) for each source rather than simply tak- 
ing the maximum likelihood estimate of the photometric redshift 
fit in constructing the redshift distribution. This has the advantage 
of giving a more robust reconstruction of the redshift distribution 
(Sheth & Rossi 2010), although it should be remembered that the 
distribution produced in this way is essentially a convolution of the 
true redshift distribution and the typical photo-z error. 

In addition to the SPIRE-mm sources, we can also estimate the 
redshifts of the 16 SPIRE-undetected mm sources, using the well- 


established radio-mm spectral index (Carilli & Yun 1999). For the 
sources in Lockman Hole North, Lindner et al. (2011) quote esti- 
mates of the redshift for all 39 1 .2-mm sources with radio identifi- 
cations, using, where possible, a combination of 1.2-mm, 1.4 GHz 
and lower frequency GMRT 325 and 610 MHz observations. We use 
their redshift estimates for the nine sources without SPIRE counter- 
parts in the Lockman Hole North field. In GOODS-N, we again use 
the spectral index approach as described by Carilli & Yun ( 1999) 
to estimate redshifts for the seven SPIRE-undetected mm sources, 
assuming the same values for a rad io and a SU bmm as Lindner et al. 
(2011), ofradio = -0.58 and a sub mni = 3.2, to ensure consistency 
between the two samples. One exception to this is G14, otherwise 
known as HDF850. 1 , for which we used the radio-far-IR photomet- 
ric redshift of z = 4. 1 ± 0.5 determined by Dunlop et al. (2004). To 
assess how reliable the mm-radio spectral index technique is for this 
sample, we compare redshift estimates produced via this method to 
the known redshifts for the 22 sources in our sample with good 
redshift information. We find a photo-z accuracy of | Az|/(1 + z) = 
0.24, or a mean difference (|Az|) = 0.6, for photometric redshifts 
estimated from the radio-mm spectral index alone. 

The full set of SPIRE-mm-radio and mm-radio photometric 
redshifts is presented in Tables 3 and 4 for the GOODS-N and 
Lockman Hole North fields, respectively. 

The redshift distribution for our SPIRE-mm sources and SPIRE- 
undetected mm sources is shown in Fig. 5. SPIRE-mm sources show 
a strong peak at z ~ 2-3, similar to other known submm/ mm sam- 
ples. Interestingly, the redshift distribution for SPIRE-undetected 
sources shows that they are predominately found at z > 3. Only 
tour of the SPIRE-detected sample are estimated to be at z > 3.5, 
while 13 of the 16 SPIRE-undetected mm sources are found above 
this redshift. This is not surprising as we have already found in 



Redshift Redshift 


Figure 5. Left: the redshift distribution of SPIRE-mm sources as determined by use of submm/mm-radio photometric redshifts. SPIRE-detected mm sources 
are shown as a grey shaded histogram. The distribution of SPIRE-mm sources with spectroscopic redshifts or reliable opticai/near-IR photometric redshifts 
from Tables 1 and 2 is shown by the dashed line. Redshift estimates from the radio-mm spectral index for sources from our parent mm sample without SPIRE 
detections are shown as a solid line. Right: the redshift distribution of SPIRE-mm sources (grey shaded) and all mm sources (grey hatched) compared to 
850 pm selected sources (solid blue line) from a compilation by Chapman et al. (in preparation), 1.1 -mm selected sources (dotted red line) from C09 and 
870 pm sources (dashed line) from Wardlow et al. (201 1). In each case, the histogram is scaled to the total number of sources in our sample. The mm-selected 
sources are seen to peak at higher redshifts than those selected at ^850 pm, with a significant excess at z > 3 made up of primarily mm-only (i.e. SPIRE 
undetected) sources. 
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Section 4.2 that the sources with the largest mm-to-SPIRE colours 
are found to be at the highest redshifts. 

Howe er, this large population of SPIRE-undetected mm sources 
at high redshift is somewhat at odds with the redshift distribu- 
tions of revious submm surveys, which are not limited by SPIRE. 
Shown for comparison is the redshift distribution of 850 pm se- 
lected sources, from a compilation by Chapman et al. (in prepara- 
tion), 1.1-mm sources in GOODS-N from C09 and 870 pm sources 
in the LESS (WeiB et al. 2009) from Wardlow et al. (201 1, hereafter 
Wll). Comparing the redshift distribution of SPIRE-mm sources 
and the submm (850/870 pm) selected samples of Chapman et al. 
(2005) and Wll, we find a clear deficit of z < 2.5 SPIRE-mm 
sources. In addition, once the mm-only sources are included, we 
find many more z > 3.5 sources than either Chapman et al. or W1 1. 
Much be:ter agreement is seen with the C09 sample, although this 
work shows somewhat fewer sources at low z (i.e. z < 1), poten- 
tially due to the smaller area probed in GOODS-N compared to our 
larger sample. A Kolmogorov-Smirnov test suggests that the C09 
redshifts and our mm-selected sample have a >99 per cent chance 
of being drawn from the same underlying redshift distribution. 

One potential reason for the disagreement between the redshift 
distribution of the mm and 850/870 pm sources is the depth of radio 
data available; in this work, we have used the two deepest blank- 
field images of the sky at 1 .4 GHz a\ ailable, with a typical rms noise 
of ~3-4 pjy. This is critical as of the 13 z > 3.5 SPIRE-undetected 
mm sources, nine are fainter than 30pJy at 1.4 GHz. C09 also 
make use of the Morrison et al. (2010) radio data in GOODS-N. In 
ECDFS, the 1.4 GHz imaging has an rms of 6.5 pJy at its deepest, 
and thus Wll will miss the counterparts for many of the highest 
redshift sources. Similarly, the Chapman et al. list is based on lower 
quality radio data over several fields, and thus is also unlikely to 
identify tne highest redshift submm sources. 

Wll attempt to account for submm sources which do not have 
robust optical/near-IR or radio counterparts by investigating the 
redshift distribution of optical-near- IR sources neighbouring their 
submm-cetected sample, and this correction is included in the red- 


shift distribution shown in Fig. 5. However, while this means that 
the W1 1 redshift distribution is not limited by the depth of the ra- 
dio (or mid-IR) imaging used for identification, this correction is 
limited by the depth of the available optical/near-IR imaging. Wll 
estimate that even after this correction, 20 per cent of LESS sources 
are still unaccounted for. 

While it appears that our work here confirms that radio-detected, 
mm-selected samples appear at higher redshift than those selected at 
^850 pm, it is difficult to determine whether this result is not simply 
a product of our use of deeper 1 .4 GHz radio data for identifications. 
To investigate the possibility that the ^-correction between submm 
(i.e. 850 pm) and mm wavelengths (i.e. 1.2 mm) is the cause for 
the difference in the redshift distribution, in Fig. 6 we plot the 
limiting IR luminosity assuming a single modified blackbody with 
dust temperature 30 and 40 K as a function of redshift. The flux 
limits at 350, 850 and 1.2 mm are taken to be 10, 5 and 2mJy, 
respectively. The 350- pm and 1 .2-mm limits were chosen to be 
representative of the sample under study here, while the 850- pm 
limit was chosen to be representative of current generation sur- 
veys at, or near, this wavelength (Mortier et al. 2005; WeiB et al. 
2009). 

From Fig. 6 it is clear that, irrespective of dust temperature, 
850 pm is marginally more sensitive in the redshift range 1 < z < 
3, while 1.2 mm is more sensitive at z > 3. We can also see that 
the 350- pm limiting luminosity increases rapidly as a function of 
redshift, crossing the 1.2-mm limit at z ~ 3. This explains the lack 
of SPIRE-mm sources seen above this redshift. 

While the difference in luminosity is small, the shape of the lumi- 
nosity function at these extreme luminosities is very steep (Chapman 
et al. 2005; Wll). At 1 < z < 2, the 850-pm and 1.2-mm sensi- 
tivities differ by ~-0.1 decades in Lj R . Using the Chapman et al. 
(2005) determination of the IR luminosity function, extrapolating 
to the limiting luminosity shown in the left-hand panel of Fig. 6 fora 
T d = 30 K SED (Lir = 12.4 Lq ; logio# ~ -5.3 Mpc -3 per decade) 
and assuming the full volume of our parent surveys between z = 1 
and 3 (0.24 deg 2 ; 5.4 x 10 6 Mpc 3 ), we could miss approximately 



Redshift Redshift 


Figure 6. Limiting IR luminosity as a function of redshift for flux-limited selection at three submm/mm wavelengths: 5350 > 10 mJy (solid line), S &50 > 5 mjy 
(dotted lins) and $ 1 . 2 mm > 2mJy (dashed line). In each case, we have assumed a modified blackbody SED with dust temperature 7 d = 30 Oeft) and 40 K 
(right). In joth panels, 850 pm is more sensitive in the redshift range 1 < z < 3, while 1.2 mm is more sensitive at z > 3. 
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five 850 pm detected sources via a 1.2-mm selection. Similarly, in 
the redshift range 3 < z < 5, assuming the sensitivities differ by 
~0.2 decades in Lir and again using the Chapman et al. (2005) 
luminosity function estimate at Lir = 12.4 Lq, we could detect an 
extra MO sources via a 1 .2-mm selection at the quoted flux limits. 
These numbers are crudely in line with the excess, and decrement, 
of sources we observe in redshift distribution of our mm-selected 
sources at 1 < z < 3 and z > 3, respectively. Performing a similar 
exercise assuming T D = 40 K would give a decrement of about one 
and excess of about two sources at 1 < z < 3 and 3 < z < 5, respec- 
tively. Thus, if the evolution of the 850 pm-1.2mm AT-correction 
is the cause of the observed differences in the redshift distribution, 
it must 3e relatively cold (T D ~30K) sources which make up the 
discrepant population. 

Of course, these calculations assume a particular optical depth 
r = (vlnuo)^, although it is worth noting that the differences in 
sensitivity between 850 pm and 1.2 mm are enhanced for lower 
values of /?, and are relatively insensitive to vq* 


5.3 Contribution of mm-selected sources to the star formation 
rate density 

Now that we have some estimate of the redshift for all of our 
parent c: radio-detected mm-selected sample, we can estimate their 
contribution to the IR luminosity density (p IR ) and hence the star 
formation rate density (SFRd) of the Universe. Given the wide 
redshift distribution of mm sources, we estimate pir in four redshift 
bins: 0.5 <z< 1.5, 1.5 <z< 2.5, 2.5 <z<3.5and3.5 <z<5.The 
redshift bins were chosen so as to match previous measurements (i.e. 
C05; W: 1) while still maintaining a reasonable number of sources 
per bin. The contribution of each source is estimated by taking 
the measured IR luminosity (Lir) and dividing by the maximum 
volume a source could be detected in (V raa x)- Incompleteness in 
the mm sample is taken into account by dividing each object by 
the completeness at its mm-wavelength flux density. Completeness 
estimates for the 1.2-mm source list in Lockman North are given 
by LincLier et al. (201 1). Unfortunately, Penner et al. (201 1) do not 
provide completeness estimates for their source list in GOODS - 
N, so v.e calculate the completeness as a function of 1.16-mm 
flux density by injecting sources at a range of flux densities into 
their combined MAMBO-AzTEC image. Sources are considered 
recovered if they are found within 3 arcsec of the injected position 
with a psak S/N of a > 3.8. 

For e^ch source, we wish to use the most reliable redshift avail- 
able. For the 22 sources with spectroscopic or good optical/near- 
IR photometric redshifts, we use these measurements, while for 
the 25 SPIRE-mm without prior redshift estimates we use our 
submm/mm-radio photo-z values from Section 5.1 . For the remain- 
ing 16 mm-ocly sources, we use estimates from the radio-mm 
spectral index, as described in Section 5.2. 

However, given the large photometric redshift errors of the 
submm/mm-radio and radio-mm spectral index estimates which 
form the bulk of our sample, we estimate p K via a Monte Carlo ap- 
proach in which realizations are produced by placing each source at 
a redshift randomly drawn from the PDF of the submm/mm-radio 
fitting. In each of the realizations, sources with good spectroscopic 
redshifts are held fixed. Sources with optical/near-IR photomet- 
ric redshifts are assumed to have Gaussian PDFs with the width 
determined by the quoted errors from Thble 2. Similarly, mm-only 
sources vith radio-mm spectral index redshifts are assumed to have 
Gaussian PDFs with width given by the typical A z found in Sec- 


Table 5. Contribution of mm-selected sources to the IR luminosity 
density. 


Redshift range 

Per 

(log io Lq Mpc -3 ) 

PSFR 

(log io Mq yr” 1 Mpc' 3 ) 

0.5 <z< 1.5 

7.1 ±0.1 

-2.7 ±0.1 

1.5 < z < 2.5 

7.72 ± 0.07 

-2.05 ± 0.07 

2.5 < z < 3.5 

7.86 ± 0.08 

-1.91 ±0.06 

3.5 < z < 5 

7.7 ±0.1 

-2.1 ±0.1 


tion 5.2 (A z = 0.6). In addition, mm-only sources arc assumed to 
be at a dust temperature of 35 K. 

In each Monte Carlo realization, we assess the contribution of 
each source to each redshift bin by taking the sum £7 4k/^L> 
where sum n runs over the sources which happen to fall in that 
redshift bin in a given realization, and L® and V ^ for source i are 
calculated using the best-fitting modified blackbody. To calculate 
Fmax> we find the minimum and maximum redshift at which an 
observed source would still be found in our sample, i.e. >3. 8a 
at 1.16 or 1.2 mm, depending on the field, and >5a at 1.4 GHz, 
taking into account variations in both the luminosity distance and 
AT -correction. The posterior probability of and the mean redshift 

in each bin are then inferred from a histogram of the Monte Carlo 
realizations, fitting a Gaussian to return the peak likelihood and 
variance. 

Based on this analysis, our best estimate of the contribution of 
mm-selected sources to the IR luminosity density is given in Table 5. 

Fig. 7 compares these values to previous measures of the contri- 
bution to the IR luminosity density by 850 pm and LESS 870 pm 
sources (Chapman et al. 2005; W1 1), and a prediction of the con- 
tribution from L® > 10 12 * 5 Lq from the model of B6thermin et al. 
(2010). At z ~ 2.5, our estimates appear lower than both the con- 
tribution from MS50 pm sources and the predictions for Lir > 
10 I2,5 Lq sources. However, there are a number of ways to ac- 
count for this. First, as we have seen in Section 5.2, mm-selected 
surveys are less sensitive than their 850 pm counterparts at 1 < z < 
3. However, mm surveys must be less sensitive in general as the 
areal density of mm-detected sources is much lower than equivalent 
850 pm samples. Including sources without radio identifications, 
our parent sample contains 82 mm-detected sources in 0.24 deg 2 , 
while the LESS survey detect 126 sources at 870 pm in a 0.25 deg 2 
and SHADES detected 120 sources at 850 pm across two fields 
totalling 0.22 deg 2 . 

Given this, it is puzzling that our estimate at z ^ 4 is somewhat 
larger than that found for 870 pm sources by Wll; however, as 
we have seen in Section 5.2, there are good reasons to expect a 
significantly larger contribution at z > 3, as mm surveys are more 
sensitive at these redshifts than their 850 pm counterparts, and we 
have the advantage of significantly deeper radio data with which to 
identify our mm sources. The Wll analysis is missing the 45 per 
cent of LESS SMGs that do not have identifications; similarly. 
Chapman et al. (2005) only have radio identifications for M0 per 
cent. Here we have radio identifications for 78 per cent of our parent 
mm-detected sample. 

It is of interest that our measurement of the SFRd at z ~ 4, while 
higher than the Wll estimate for 870pm sources, is still some- 
what short of that predicted for IR luminous galaxies (logioL^ > 
12.5 Lq) by B6thermin et al. (2010), and MOO lower than the 
total SFRd estimated from a combination of UV, radio and IR 
measurements (Hopkins & Beacom 2006). Thus, while select- 
ing sources at mm wavelengths offers an effective way to isolate 
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Figure 7. Contribution of mm-selected sources to the 1R luminosity density 
evolution. For comparison, we show the estimates of Chapman et al. (2005) 
and Wll at z = 0-4. As ECDFS is underdense at submm wavelengths, 
W1 1 present two estimates: one based on the observed number of sources in 
ECDFS ai d another which has been scaled to the SHADES 850 pm num- 
ber density (Coppin et al. 2006). Also shown is the predicted contribution 
by > 10 12 * 2 Lq ultraluminous IR galaxies to the IR luminosity density from 
the empirical model of B&hermin et al. (2010). At z < 3, our estimate is 
significantly lower than those found in both previous studies and models, 
while hig/:er than previous studies at z ~~ 4. These disagreements can be 
explained when considering the flux limits and biases of the comparison 
samples. 


luminous starbursts at high redshift ( z > 3), it is clear that mm- 
selected sources, at least at the imaging depth possible with current 
facilities, do not make a significant contribution to the SFRd at any 
redshift. 

Due to the strong dust temperature bias, mm-selected sources 
never cor tain the bulk of IR luminous sources at any redshift. Taking 
the Bdthermin et al. (2010) prediction as a guide, our mm-selected 
sample nsver makes up more than ^50 per cent of the expected 
SFRd of IR luminous sources (logi 0 LiR > 12.5). We can see the 
reason fcr this clearly in Fig. 6: at T D = 30 K we could expect 
our mm sample to reach this level, while at T D = 40 K we are only 
sensitive to IR sources above a luminosity of log ]0 > 12.7. Thus, 
for the Bethermin et al. (2010) prediction to hold, this shortfall must 
be made up of IR sources with logio Ur < 12.7 and T D > 30 K. 


6 CONCLUSIONS 

We have presented an analysis of 63 radio-detected, mm-selected 
sources. Our parent sample of mm sources was drawn from the re- 
cent MAMBO 1 .2-mm imaging of Lockman Hole North by Lindner 
et al. (201 1) and the combined MAMBO- AzTEC 1.16-mm image 
of GOODS-N presented by Penner et al. (2011). Radio identifi- 
cations a.e found for the Penner et al. (201 1) sample, with coun- 
terparts for 24 of the 41 1.16-mm sources presented. These new 
identifications were combined with the 39 (out of 41) radio iden- 
tifications presented by Lindner et al. (2011) for their MAMBO 


sample to give a sample of 63 radio-detected mm sources. We have 
cross-matched this mm-selected sample with SPIRE data from the 
HerMES project, making use of a map-based technique which re- 
duces the problem of confusion in the SPIRE data by using radio 
and 24 ^im positions as a prior. The main results of this work are as 
follows. 

(i) After comparing with the SPIRE data, we found 47 mm 
sources with a detection in one or more SPIRE band, giving a 
detection rate of 74 per cent. 

(ii) For SPIRE-detected mm sources with spectroscopic or good 
quality optical/near- IR photometric redshifts, we found a tight cor- 
relation between dust temperature (T D ) and IR luminosity (L K ). 
This correlation is found to be offset to colder dust temperatures 
from the local relationship (e.g. Chapman et al. 2003). 

(iii) Using the observed tight relationship between dust temper- 
ature (T d ) and IR luminosity (Lir) to break the degeneracy between 
redshift and dust temperature, we have demonstrated that photo- 
metric redshifts from submm/mm data alone offer an accuracy of 
| Az|/(1 4- z) = 0.16 ((| Az|) = 0.51). If constraints from the radio- 
far-IR correlation are included, this is improved to | Az|/(1 -f - z) = 
0.15 «|Az|) = 0.45). 

(iv) Via a combination of spectroscopic, optical/near- IR pho- 
tometric and submm/mm-radio photometric redshifts, the redshift 
distribution of mm-selected sources was presented. We found that 
the redshift distribution of mm-selected sources peaks at higher 
redshifts than those selected at ^850 |im, with a long tail to high 
z (z > 3.5) made up of primarily mm-only (i.e. SPIRE undetected) 
sources. This difference may be explained by a combination of the 
depth of the radio data used to identify the sources, as well as the 
evolving /^-correction between submm and mm wavelengths. Good 
agreement was found between this study and previous mm-selected 
samples. 

(v) We measured the contribution of mm-selected sources to the 
SFRd of the Universe. We found that mm-selected sources do not 
contribute significantly to the global SFRd, nor do they make up the 
bulk of IR luminous sources, at any redshift. 
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